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Abstract 



The Office of Naval Research Tropical Cyclone Motion 
initiative is a five-year program to improve basic 
understanding of tropical cyclone motion. On 27-30 January 
1988, a meeting was held in Monterey, California to: (i) 

Review the first year's research activities; (ii) Discuss 
issues and plan future research; and (iii) Begin discussion 
of tentative hypotheses that might be explored in a field 
experiment in the western North Pacific region during summer 
1990. Each of these topics will be briefly discussed to 
indicate the progress and plans of the initiative. 
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1 . 



Introduction 



A five-year basic research program to improve 
understanding of tropical cyclone motion began 1 October 
1986 under the sponsorship of the Office of Naval Research 
Marine Meteorology Program (R. F. Abbey, Jr, Program 
Manager) . This program involves theoretical studies, 
analysis of existing observational data, and a field 
experiment in the western North Pacific region during summer 
1990. A series of workshop reports (Elsberry, 1986; 1987a; 
1987b) describe respectively; the planning of theoretical 
studies; possible observing systems for tropical cyclone 
studies; and a reassessment of the program in view of 
elimination of aircraft reconnaissance in the western North 
Pacific during 1987. 

A workshop was held on 27-30 January 1988 to review the 
progress in the first year of the initiative. A list of 
attendees is given in Appendix A. We continue to benefit by 
the participation of cooperating agencies, such as the 
Hurricane Research Division (HRD) of the National Oceanic 
and Atmospheric Administration (NOAA) , the U.S. Air Force, 
and the Naval Environmental Prediction Research Facility. 

The agenda for the workshop is provided in Appendix B. 
The objectives of the workshop were three-fold: (i) Present 

reports of progress during the first year; (ii) Discuss 
the status of the issues being addressed and identify future 
areas of research; and ( iii) Discuss potential hypotheses 
for the 1990 field experiment. 
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R. Abbey of ONR opened the meeting by briefly 
describing the list of current contractors on the project 
(Table 1) . Analytical studies are defined here to include 
observational analyses and primarily numerical approaches 
are grouped with theoretical studies. Although each 
contractor is listed in only one category, most 
investigators actually contribute in more than one category. 

The first day and part of the second day of the 
workshop were devoted to progress reports by ONR contractors 
and by participants from cooperating agencies (see agenda in 
Appendix B) . Although some of this work is nearing 
publication, other portions are recent results that may 
change with further investigation and analysis. 

Consequently, no attempt is made to describe these 
presentations in detail. Some of the presentors have 
provided summaries that are included in Appendix C. 



Table 1. Contractors participating in the ONR Tropical 
Cyclone Motion Accelerated Research Initiative (provided by 
R. F. Abbey, Jr.) 

Experimental 

V. Lally, National Center for Atmospheric Research 
Analytical Studies 

R. Elsberry, Naval Postgraduate School 

W. Frank, Pennsylvania State University 
W. Gray, Colorado State University 

G. Holland, Bureau of Meteorology Research Center 

R. Merrill, University of Wisconsin 
C. Ramage, University of Hawaii 

T. Schroeder, University of Hawaii 

T. Tsui, Naval Environmental Prediction Research Facility 
Theoretical Studies 

S. Chang, Naval Research Laboratory 

J. Chi, University of the District of Columbia 

T. Krishnamurti , Florida State University 
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B. Morton, Monash University 
W. Schubert, Colorado State University 
R. Smith, Monash University 
B. Wan, University of Hawaii 
T. Williams, Naval Postgraduate School 

2. Discussion Sessions 

Three sessions on the second day were devoted to open 
discussions of issues related generally to theoretical 
studies (R. T. Williams, discussion leader) , observational 
studies (G. J. Holland, discussion leader) and numerical 
studies (W. Schubert, discussion leader) . Since the format 
was left open, these discussion leaders used a combination 
of personal research with results presented in the progress 
reports to highlight achievements and raise questions/issues 
for discussion. Rapporteurs had been assigned to summarize 
these discussions. The summaries prepared by these 
rapporteurs are included below. 

a. Theoretical Studies [H. Willoughby, Rapporteur] 

(1) Accomplishments : Theoretical investigations under 

this initiative have focused on motion of barotropic 
nondivergent and divergent vortices. On an f plane, a 
spatially uniform environmental current simply advects the 
vortex with no distortion, provided that it induces little 
divergence. A rotating imposed divergence can induce 
trochoidal motions if it moves with the swirling wind where 
it is imposed. On a beta plane without environmental flow, 
the vortices' outer part plays a greater role than their 
cores in determination of the northwestward or westward 
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motion that arises. Vorticity advection — expressed as a 
balance among radial advection of axisymmetric vorticity by 
the asymmetries, tangential advection of the asymmetric 
vorticity by the mean swirling flow, and meridional 
advection of planetary vorticity also by the mean flow (beta 
effect) — is the primary physical process responsible for 
the motion (Fig. 1) . This balance suggests that the 
asymmetries are Rossby waves that propagate on the radial 
gradient of relative vorticity in the cyclone and that the 
beta advection term, which is smaller than the relative 
vorticity advection, acts to force the waves. Although the 
linear theory of Willoughby (1987) seems to describe the 
westward drift and pattern of asymmetry of a vortex on a 
beta plane, an unreasonably fast motion to the northwest is 
predicted, which indicates that nonlinearity is important. 

The flow in a moving vortex may be partitioned into the 
axisymmetric vortex, the relative environmental current, and 
the Rossby- or gravity-wave asymmetries. It is the 
wavenumber one asymmetry that influences the vortex motion 
most strongly. A horizontally shearing environmental 
current excites wavenumber two, which may interact 
nonlinearly with wavenumber one to influence both wavenumber 
one and the motion. Gradients of potential vorticity, which 
may be due either to horizontal variations of the 
environmental vorticity or to topography of the free surface 
that sustains the environmental flow, may affect the motion 
much as the beta effect does. More complicated and 
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Fig. 1 Streamfunction tendency (m 2 s“ 2 ) after 12 h in the 
nondivergent barotropic model simulations of Fiorino (1987) 
due to the (a) linear 0v term; (b) advection of symmetric 
vorticity by the asymmetric flow; and (c) advection of 
asymmetric vorticity by the symmetric flow. Positive 
(solid) and negative (dashed) tendencies correspond to 
increasing anticyclonic and cyclonic circulations 
respectively. Since the contour intervals are 0.02, 0.10 
and 0.02 respectively, the leading term is the advection of 
symmetric vorticity by the asymmetric flow, which is 
illustrated in panel d by the streamlines and relative 
vorticity contours (solid) after 72 h. Notice that the 
large-scale anticyclonic (cyclonic) gyres to the right 
(left) of the central region of high vorticity establish the 
asymmetric flow. 
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realistic environments would include subtropical 
anticyclones. Tropical Upper Tropospheric Trough (TUTT) 
cells, extratropical troughs and other tropical cyclones. 

Vortex asymmetries, which are generally a pair of 
counterrotating gyres, occur in two kinds: the beta gyres, 

which are induced by advection of planetary vorticity or 
other synoptic-scale influences and are comparable in 
horizontal scale with the vortex; and the alpha gyres, which 
are confined to the vortex center and may be eliminated by a 
suitable redefinition of coordinates. Presence of alpha 
gyres makes the total-flow streamlines take the form of 
eccentric circles about the pressure center. In a 
baroclinic model, the alpha gyres may reflect a slope of the 
rotation axis. Minimization of the Lagrangian, which 
entails minimization of the amplitude of the alpha gyres f 
provides a basis for determination of vortex motion in a 
linear model. 

(2) Opportunities for future research : Nonlinear 

processes — interactions among asymmetric features and 
between the asymmetries and the axisymmetric vortex — must 
be important to tropical cyclone motion, as must matching of 
the relative-vorticity Rossby waves in the vortex core onto 
planetary Rossby waves in the surroundings. The 
interactions and matching may involve inward communication 
of environmental influences or outward radiation of energy 
and angular momentum, perhaps acting to dissipate the mean 
vortex. In addition to nonlinear processes, linear 
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barotropic damping of the asymmetries by the horizontal 
shear in the mean flow may occur and asymmetric surface 
friction due to the vortex motion may excite wavenumber-one 
asymmetries that feed back on the motion. 

Rossby-wave asymmetries in a stratified barotropic 
model without vertical shear have reduced equivalent depth 
and Rossby deformation radius, so that the beta effect 
should be reduced. With vertical shear, differential 
advection of the axisymmetric flow by the relative 
environmental flow and twisting of the horizontal 
environmental vorticity by the axisymmetric secondary 
circulation should induce asymmetries, and perhaps lead to 
inclination of the mean rotation axis. 

Future research should focus on solution matching at 
the vortex periphery, frictional effects and baroclinity. 
b. Observational Studies [J. C.-L. Chan, Rapporteur] 

(1) General : The observational studies in this 

initiative have produced only some preliminary results. 

These results are significant for two reasons. First, they 
confirm some of the theoretical or numerical results. 

Second, some important discrepancies exist between these 
observational findings and those from the theoretical and 
numerical investigations. In this brief summary, these 
agreements and discrepancies are described. Questions 
raised during the discussion session are also summarized. 
Finally, future work to address these questions is proposed. 
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(2) Structure of the beta-gyres : Numerical and 

theoretical studies have shown that nonlinearity in a 
barotropic framework leads to the formation of asymmetric 
gyres on opposite sides of the vortex (Chan and Williams, 
1987; Fiorino, 1987; Willoughby, 1987; Williams and Peng, 
workshop presentation) . The gyres at radii 500-800 km from 
the center of the vortex that arise from the beta-effect are 
referred to as the beta-gyres. Observational studies (Chan, 
1986; Holland, workshop presentation) have documented the 
existence of the beta-gyres that have the same spatial scale 
as those predicted in the theoretical studies. Holland 
found that these gyres appear to be oriented east-west 
across the vortex, develop after the cyclone forms and 
disappear as the cyclone decays. 

The major issues that need to be addressed here include 
the structure and orientation of these gyres. How does the 
structure evolve with time or vary with different 
environmental and storm-related parameters such as size and 
intensity? Is the orientation of these gyres necessarily 
east-west? In Chan's (1986) study of Supertyphoon Abby, the 
gyres seemed to rotate as the heading of Abby changed from 
north-northwest to north-northeast. Fiorino (1987) also 
found the gyres in his model rotate as the vortex moved 
first toward the north and then toward the northwest. 

Another question is whether the "gyres” identified in 
the observational studies are simply part of the large-scale 
environment. Holland's observations that these gyres seem 
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to develop and move with the storm suggests that they are 
not simply due to the environment. One could resolve this 
issue if the large-scale environment could be identified 
clearly. However, the consensus during the workshop appears 
to be that a unambiguous definition of the "large-scale" 
remains elusive. 

To substantiate further the existence and document the 
structure of these gyres, data from storms of different 
sizes and intensities (including monsoon depressions) need 
to be analyzed. The data sources include the Australian 
Monsoon Experiment (AMEX) data set and the reconnaissance 
aircraft data processed by Colorado State University. In 
the latter set, both composite and individual case studies 
should be attempted. The geographical as well as the storm- 
direction oriented coordinate systems should be used, but 
with the symmetric component of the vortex subtracted. The 
data can also be stratified based on synoptic flow patterns, 
such as the subtropical ridge or even derived quantities 
such as relative vorticity or its gradient. A spectral 
decomposition of the azimuthal eddy would provide the 
amplitude and phase of the wavenumber one asymmetry. 

(3) The westward drift : In the numerical study of 

Fiorino (1987), a discrepancy exists between the vector of 
the vortex movement and that of the "ventilation flow" 
(asymmetric flow averaged within 300 km of vortex center) . 

In his presentation at the workshop, L. Carr of the Naval 
Postgraduate School showed that the difference between the 
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two vectors consistently is about 0.3 m/s towards the west. 
In the linear model of Chan and Williams (1987) , the 
displacement of the streamfunction minimum from the relative 
vorticity maximum is also about this order of magnitude. 

Carr also examined observational studies of George and 
Gray (1976) , Chan and Gray (1982) and Holland (1984) in 
which storm translation was compared with a "steering flow", 
which was represented by the average around an annulus at a 
specific radii from the storm center. The difference 
between these two vectors is generally westward with a 
magnitude of about 1-2 m/s. This difference vector also 
rotates slightly and varies in magnitude for different 
stratifications. Holland reported that the difference 
vector is about 1-2 m/s in the AMEX data set, compared to 
the 0.3 m/s magnitude in the numerical experiments of 
Fiorino (1987) . 

The question is why a discrepancy exists between the 
numerical model and the observational results. Is it 
because of a data problem — perhaps the annulus over which 
the flow is analyzed was not chosen correctly? Or is it 
because the model is too simple? What is the magnitude of 
this difference vector in other data sets? 

Thus, the steering needs to be recomputed using 
different annuli and for different levels and deep-layer 
means with all available data sets. The workshop 
participants agreed that at least the flow at 850, 700, 500 
and 200 mb as well as the 850-300 mb deep-layer mean should 
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be analyzed. New data sets should include rawindsondes 
composites, reconnaissance aircraft data (both individual 
case and composites) , AMEX data set and the Hurricane 
Research Division (HRD) data (including dropwindsondes) . 
Composites can also be made based on vorticity gradient 
patterns to determine whether the discrepancies arise as a 
result of environmental vorticity gradients. 

"Steering" flow also can be computed from numerical 
models using the same definition as in the observational 
studies. This can be done in the numerical experiments with 
and without basic flow included. 

(4) Vortex profiles : In many observational studies 

(composites and individual cases) , the tangential winds 
remain cyclonic out to 15° latitude radius unless the upper- 
level flow is included. Thus, observed tropical cyclones do 
not have a tendency to become angular momentum neutral as 
was suggested by Willoughby (1987) and Fiorino (1987). 
Further, although the relative vorticity becomes 
anticyclonic away from the storm center, it does so at radii 
that are at least a few times the radius of maximum winds. 

On the other hand, in most of the vortices used in numerical 
studies of DeMaria (1985), Chan and Williams (1987) and 
Fiorino (1987) , the relative vorticity becomes anticyclonic 
at twice the radius of maximum wind. 

The major question is therefore whether the numerical 
results are an outcome of the "unrealistic" specification of 
the vortex profile. What effects will the anticyclonic 
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relative vorticity produce in the movement of the vortex? 

As observed vortices have positive tangential winds at such 
a large radius, is the concept of a "bounded” vortex 
realistic? 

To address these questions, the models should be 
initialized with more realistic vortices. Vortex profiles 
in different basins can also be compared to determine their 
inherent differences (e.g., western Pacific tropical 
cyclones typically form in a monsoon trough whereas Atlantic 
cyclones form in a trade-wind regime) . 

(5) Summary : In addition to the three main topics 

discussed in the session, some potential new observing 
systems (e. g. , satellite microwave sensors, Doppler radars, 
drifting buoys, etc.) that might be provided by the Air 
Force or the Navy were also described by C. Guard. They 
point to a significant increase in the amount of available 
data if these systems become operational by the time of the 
field experiment. 

It appears that much more can be done with the existing 
data sets to either verify or point out the deficiencies of 
the model results. They also can provide some better 
insights into some of the physical processes that are 
responsible for tropical cyclone motion, 
c. Numerical Modeling Studies [M. Fiorino, Rapporteur] 

(1) Accomplishments : The major accomplishments during 

the first year of the initiative were: (i) An understanding 

of the role of the large scale flow and vortex structure in 
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tropical cyclone motion using nondivergent and divergent 
barotropic prognostic and diagnostic models (Fiorino, 
Williams, Willoughby, Chan and DeMaria) ; (ii) Upper-level 
flow influences using a limited-area primitive equation 
model (Chang) ; (iii) Real-data forecasts using a spectral 
global model (Krishnamurti) ? and (iv) Predictability of 
tropical cyclone motion (DeMaria) . 

Following a discussion of these accomplishments, Dr. 
Krishnamurti was asked to review the work at Florida State 
with global models and to provide additional results. He 
commented on the predictability limits derived by DeMaria 
based on a presumed observation spacing of 3000 km in the 
Atlantic. While this area has virtually no upper air 
observations east of 60° W to the African coast, the 
combination of many aircraft reports, satellite winds in the 
lower and upper levels, and ship reports, combined with a 
good 4-D assimilation system, tends to reduce the 
"effective" station separation. Thus, there may be optimism 
for accurate model predictions despite the sparsity of 
rawindsonde observations. Higher resolution may be needed 
in the Atlantic compared to the Pacific, although 
climatological or environmental reasons for this apparent 
requirement are not known. Krishnamurti also discussed the 
development of a "reverse" Kuo scheme in which satellite 
rainfall estimates are used to recover a 
thermodynamic/dynamic structure that will support the 
observed rainfall. The intent is to force the model to 
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produce typical tropical rainfall rates at the beginning of 
the integration rather than relying on the model to generate 
rain from a nearly nondivergent initial flow. This spinup 
usually takes 18-36 h to build the precipitation to 
reasonable levels. 

J. Chan commented on the performance of the European 
Center for Medium-range Weather Forecasts (ECMWF) model for 
typhoon Wayne of 1986. This cyclone made four loops in the 
South China Sea and the ECMWF model predicted two of the 
loops. In these two cases, the large-scale forcing appeared 
to be stronger. In the two cases in which the model failed 
to forecast a loop, the storm was weak or the synoptic flow 
weak. This suggested a vortex-structure effect that could 
not be handled due to the inability of the analysis to 
resolve cyclone scales. 

M. DeMaria will reanalyze the cases in his 
predictability experiments in which there were large 
sensitivities to the initial conditions and/or synoptic 
flow. 

(2) Future questions : 

(i) Can theoretical studies of cyclone-generated 
asymmetries provide initial conditions for numerical 
forecast models? 

(ii) What numerical experiments can be used to 
quantify the importance of baroclinic effects, friction, 
convection, etc. , on tropical cyclone motion? 



16 



(iii) What model initialization and data assimilation 
techniques can be used to reconstruct the poorly observed 
vortex scales (e.g., the Adjoint method)? 

(iv) How can models be designed to use wind profiler 

data? 

These questions stimulated a discussion on the future 
role of numerical modelling in understanding tropical 
cyclone motion. Real-data global simulations provide an 
opportunity to test motion concepts developed with simpler 
numerical models and analytical treatments. In this 
scenario, the global model output represents a 4-D "data 
source" for comparisons with further work and rawindsonde or 
aircraft data sets. 
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3. 



Tentative Hypotheses 



a. General Recurvature Forecast Problems [Provided by G. 

Holland] 

Tropical cyclones in potential recurvature situations 
are some of the most difficult to forecast and these 
situations produce a significant percentage of the largest 
forecast errors (for example, Sandgathe, 1987; JTWC Annual 
Typhoon Reports; Bureau of Meteorology, 1978). A few 
schematic examples in the western North Pacific region from 
Sandgathe (1987) are shown in Fig. 2. Some systems 
unexpectedly turn into the subtropical ridge. Sometimes the 
ridge builds westward ahead of a recurving cyclone, which 
instead continues on a westward track. Cyclones approaching 
a midlatitude trough have been observed to recurve, to 
continue westward, or to 'step' or loop before returning to 
a westward movement. 

These examples indicate the importance of understanding 
the interactions between cyclones and both the subtropical 
ridge and the environmental (midlatitude or TUTT) troughs. 
Unfortunately, much of our present knowledge consists of 
anecdotal accounts and empirical rules of thumb. In the 
following sections, working hypotheses are proposed. A 
structured research and field program aimed at improving our 
knowledge of these interactions is required (see Section 4) . 
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(C) EXTEND 




(d) RECURVE 




Fig. 2 Tropical cyclone motion during cyclone - subtropical 
ridge interaction for (a) normal motion, (b) through ridge, 
and (c) extending ridge, and during cyclone - midlatitude 
trough interaction for (d) recurvature, (e) stepping and (f) 
looping (after Sandgathe, 1987) . 
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(1) Tropical Cvclone/Subtropical Ridge 
Interactions : 

The aim here is to investigate the interactions that 
occur between the cyclone and the subtropical ridge that can 
lead to modification of one or both systems and thus affect 
the cyclone motion. Two specific working hypotheses are 
proposed: 

(a) The beta gyres arising from Rossby mode dispersion 
of the cyclone in a gradient of earth vorticity can 
reach sufficient amplitude to modify the subtropical 
ridge. A schematic is shown in Fig. 3. 

The issues to be addressed include the conditions in which 

significant amplitude gyres are generated, the conditions 

that affect their orientation, and the track modifications 

that result. 

(b) Thermal and indirect circulation effects 
associated with the tropical cyclone outflow jets can 
produce a westward extension of the mid- to lower-level 
subtropical ridge and cause a continued westward motion 
of the cyclone. 

The requirement here is to determine the conditions that 
generate and orient the outflow jets and the manner in which 
these jets can affect the lower troposphere. 

(2) Tropical Cvclone/Environmental Trough Interaction : 
The focus here is to investigate the interactions that 

occur when a tropical cyclone comes into close proximity to 
a progressive midlatitude trough in the westerlies or a TUTT 
cell, and to determine the types of interactions that are 
associated with the track types shown in Fig. 2. The 
specific working hypothesis is: 
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Fig. 3 Schematic of the beta gyres (dashed) as in the 
isolated vortex model of Fiorino (1987) overlaid on a 
subtropical ridge (STR) . The anticyclone (A) - cyclone (C) 
couplet can potentially move the cyclone toward the ridge 
axis and also modify the ridge (provided by G. Holland) . 
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The approach of an environmental trough may excite 
barotropically unstable azimuthal waves in the 
tropical cyclone outflow layer that extend downward 
into the lower troposphere and affect the cyclone 
motion. A schematic is shown in Fig. 4. 

The AMEX analyses by Holland show evidence of barotropically 

unstable growth that causes an equatorward extension of the 

mid-latitude trough and affects recurvature. Additional 

studies are needed to determine the conditions under which 

barotropic instability is generated and released, the types 

of perturbations that result, and the manner in which they 

can affect the cyclone motion. 

b. General Track Turning Problems [Provided by R. L. 
Elsberry ] 

A related topic to the recurvature forecast problem 
discussed above by G. Holland is the general topic of "how 
tropical cyclone turning motion is related to a time- 
changing environmental flow". As indicated in Fig. 2, a 
long-lasting tropical cyclone in the western North Pacific 
may have multiple interactions with midlatitude troughs 
prior to the actual recurvature. This leads to a "stepping" 
track in which the tropical cyclone "feints" a turn toward 
the north and then returns to a westward track. 

Considerable literature exists regarding environmental 
influences, but conventional or inner-core reconnaissance 
observations do not appear to be adequate to understand (or 
forecast) the nonlinear processes that are likely to be 
involved. It is possible that the capability to remotely 
sense or to directly observe (via a dropwindsonde survey) 
the surrounding environment can now be exploited to improve 
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Fig. 4 Cyclonic midlatitude trough impinging on the 
anticyclonic outflow of the tropical cyclone may change the 
sign of the vorticity gradient and lead to excitation of 
barotropically unstable developments shown by the dashed 
line extending equatorward (provided by G. Holland) . 
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understanding. The questions are: Where are the critical 

areas that need to be observed? What observations or 

diagnostic calculations need to be made to distinguish 

unambiguously these turning motions? 

Fiorino (1987) indicates that tropical cyclone turning 

motion (associated with •’external" (300 - 800 km) structure 

changes) is preceded (accompanied?) by reorientation of the 

spatially-averaged current. Sherman (1988) has examined 

gradients of relative vorticity (derived from operationally 

analyzed fields) in the environment of "stepping tracks." 

He finds that the tropical cyclone track tend to be oriented 

along (versus across) "channels" or "primary nodes" in the 

relative vorticity contours surrounding the tropical 

cyclone. These patterns may be related to the beta-gyres as 

in Fig. 3 or other interactions between the tropical cyclone 

and the surrounding flow. The operationally-analyzed fields 

are quite noisy and the delay in the time until the tropical 

cyclone turns is uncertain at this stage. Nevertheless, a 

tentative hypothesis might be stated: 

Tropical cyclone turning motion, or the absence 
of any turns in the immediate future, may be 
monitored by accurate representations of the 
large-scale vorticity patterns in the environ- 
ment of the cyclone. 

Additional studies are needed to determine whether 
Sherman's preliminary relationships between the vorticity 
features and storm turning are maintained. However, the 
coarse resolution and questionable data coverage in the 
tropics are not expected to allow conclusive tests. Some 
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higher resolution analyses may be available from a more 
advanced numerical model being developed at NEPRF. In 
addition, the time-dependent " forcing" fields can be 
extracted from the history tape of the model predictions of 
recurvature and non-recurvature tracks. Even if these 
studies based on analyses and predictions seem to 
substantiate the hypothesis, a field experiment will be 
required to obtain the specific observations that 
demonstrate conclusively the physical relationships, 
c. Interactions with TUTT Cells [Provided by C. Guard] 

The western North Pacific basin has many special 
characteristics that affect the development, intensity and 
movement of tropical cyclones. For example, monsoon surge 
interactions may be related to genesis, to intensification, 
and to downstream influences on motion. 

Another special influence (discussed in part by G. 
Holland above) may occur when intense TUTT cells extend 
downward into the mid-troposphere. Satellite imagery 
suggests that very active monsoon regimes also may be 
triggered by a sequence of TUTT cells aloft. A number of 
possible track effects are possible due to the interactions 
with the TUTT cells. Sadler (1978) has described multiple 
storm developments in active trough regimes. The intensity 
changes in association with outflow jets interacting with 
TUTT cells have been studied by Sadler, Merrill, Holland and 
others. TUTT cells that are within 10-12° latitude west or 
poleward of a tropical cyclone may induce rapid deepening. 
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This deepening may alter both the intensity and outer 

strength of the cyclone and the outflow jet intensity, which 

may change the environmental flow around the tropical 

cyclone. Thus, interactions with TUTT cells may lead to 

track deviations. A preliminary hypothesis is: 

Significant tropical cyclone track direction changes 
(such as steps in Fig. 2) can occur during periods 
of direct or indirect (e.g., induced monsoon surges) 
interaction with TUTT cells. 

Additional research is necessary to document potential 
track steps as the TUTT cells move past. Questions such as 
the amount of track deviation as a function of TUTT cell 
intensity, closest point-of -approach, and critical heights 
gradients need to be resolved. 

d. Relation of Model/ Analysis Resolution to Track Error 
[Provided by T. Krishnamurti] 

Very high resolution (about 50 km) global model 
simulations suggest that tropical cyclone motion can now be 
addressed in a real-data environment. The summer 1990 field 
experiment may provide adequate observations to test the 
effects of numerical model and analysis resolution on the 
accuracy of the tropical cyclone prediction. Given adequate 
observations, the track prediction accuracy should improve 
with increased numerical model resolution. Specifically, 
the hypothesis is: 

Given adequate four-dimensional observations to 
define the initial conditions, the 24-h track 
prediction error should not exceed five grid 
points of the model resolution. 

This hypothesis might be testable if sufficiently high 

resolution observations are available to define both the 
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inner-core region and the environment in which the cyclone 
is embedded. 

e. Cyclone Turning and Recurvature [Provided by W. Gray] 

A number of studies have been carried out with 
rawindsondes composites of turning storms or recurving 
tropical cyclones. New studies are in progress with the 
aircraft reconnaissance data from the western North Pacific 
These studies have not proceeded to the stage where a 
specific hypothesis can be stated. Nevertheless, it does 
appear that distinctive large-scale changes occur in these 
fields (and perhaps also in the satellite data) in advance 
of the cyclone recurvature or turning motion. In the 
future, a hypothesized sequence of events that accompany 
recurvature will be proposed for test and evaluation in the 
field experiment. 

4. Preliminary Field Experiment Design Considerations 

[Provided by G. Holland] 

The above hypotheses have been developed from recent 
work in the ONR Tropical Cyclone Motion research initiative 
Further research is continuing using numerical models, the 
AMEX data, and case studies of recurving cyclones in the 
western North Pacific Ocean. This research can be expected 
to lead to further refinement of the hypotheses and perhaps 
development of new hypotheses. 

An integral component of this research must be a 
testing of the final hypotheses with real data. Because of 
the highly transient and localized nature of many of the 
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interactions described above, routine data collection will 
not be sufficient, and a special field experiment is 
required. The western North Pacific Ocean during August and 
September is a logical site for this experiment for the 
following reasons: (i) this is the major region of interest 

to the ONR program and has a substantially different 
environment relative to other cyclone regions; (ii) the 
frequent cyclone activity there provides a high probability 
of successful data collection in a number of cyclones; and 
(iii) the current observing infrastructure and available 
additional sites provide an excellent basis for an 
experiment. 

It is proposed that the experiment be based on a fixed 
network of upper-air sounding sites that consists of an 
enhanced experimental array embedded within the current 
operational network. Additional support would be from 
mobile aircraft-based platforms, from satellite, from 
ground-based radar, and perhaps drifting buoys, 
a. Upper-air Observing Network 

The proposed upper-air observing network (Fig. 5) 
consists of a regular array of experimental stations 
embedded within the routine station network. It is 
envisaged that the experimental sites will consist of a 
combination of upgraded operational stations, portable 
sounding stations using either Loran or Omega, and vertical 
wind profilers. 
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Fig. 5 Example of a possible experimental upper-air network 
(large circles) superposed on the present observing network 
according to Sandgathe, 1987. 
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The experimental network in Fig. 5 has been designed to 
be as regular as possible, to take advantage of the high 
density sites over Asia, and to cover the main regions in 
which the interactions described above will occur. For 
example, the meridional station chains can provide excellent 
location and structural descriptions of the subtropical 
ridge, and of the very large scale structure of the tropical 
cyclones. Further, good observations would be available for 
subtropical trough/cyclone outflow interactions in the 
vicinity of Japan. The regular nature of the network 
enhances the capacity for analysis on a rectangular grid 
with a minimal loss of information. Subsets of the network 
also could be used for detailed diagnostic calculations 
using line integral or other techniques. 

This fixed network provides the basis for the 
experiment and should operate during all cyclone events for 
at least a two month period. It should be supported by ad 
hoc data collected from ships and aircraft of opportunity, 
and by sounding and remote sensing data from polar orbiting 
satellites and by the Japanese Geostationary Meteorological 
Satellite. The value of these data would be enhanced 
considerably if mobile platforms also were used to collect 
very high density data in specific regions of interest, 
b. Mobile Platforms 

Two types of mobile platforms are desirable. The first 
is aircraft with Doppler or dropwindsonde capacity to locate 
the tropical cyclone and monitor the structural changes 
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within a few hundred kilometers of the center. Although the 
cyclone intensity does not seem to be important for motion, 
there is evidence that the structure, including patterns of 
precipitation on scales of 100-130 km, can be very 
important. Further, these aircraft could be utilized to 
monitor the detailed structure of (and changes in) the 
subtropical ridge during interactions with tropical 
cyclones. 

The second mobile requirement is for high flying 
aircraft with dropwindsonde and remote sensing capacity that 
can monitor the outflow layer wind structure and 
interactions with subtropical troughs, with TUTT cells and 
with the subtropical ridge. These aircraft would need to be 
capable both of taking direct observations within the 
outflow regimes and providing high resolution dropwindsonde 
coverage from the upper troposphere. 

5. Conclusions 

The intent of this brief report has been to indicate 
the progress during the first year of the ONR Tropical 
Cyclone Motion initiative. Regular journal articles are 
being submitted on a number of topics, and new research 
areas are developing continually. It was particularly 
pleasing to see the confluence of thought on theoretical 
approaches, and the cross-fertilization with the 
observational ists based on the theoretical results. 

Finally, some tentative hypotheses for the 1990 field 
experiment are arising. Much remains to be accomplished to 
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put these in the form of testable hypotheses, and to design 
the appropriate field experiment within the severe financial 
constraints. 

The next ONR workshop is planned for 29 June - 1 July 
1988 in conjunction with the tropical meteorology conference 
in Brisbane, Australia. Scientists attending that 
conference are encouraged to attend the ONR workshop. In 
addition to presenting progress reports, refinement of 
hypotheses will be a key objective. This will contribute to 
the preparation of an operations plan, which will be 
published in 1989. 
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FIRST-YEAR REVIEW 

27-30 January 1988 
Monterey, California 

AGENDA 



27 January 1988 
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ONR Status Report - R. F. Abbey, Jr. 
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M. Fiorino - Effect of vortex structure on 
tropical cyclone motion 
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S. Chang - Model -generated structure of tropical 
cyclones 

M. DeMaria - Predictability of tropical cyclone 
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Discussion Leader; H. W. Willoughby, Rapporteur) 

- Status of the theoretical studies 
Lunch 



36 



1245 



1445 

1500 



1630 

Friday, 

0830 

0900 

0930 

1000 

1015 

1045 



1145 

1245 

1500 

1515 

1630 

Saturday. 

0900 

1030 



Discussion session #2: (G. J. Holland, Discussion 

Leader; J. C.-L. Chan, Rapporteur) - Status of 
Observational Studies 
Break 

Discussion session #3; (W. Schubert, Discussion 

Leader; R. Hodur, Rapporteur) - Status of 
Numerical Prediction Model Studies 
Close 

January 1988 

Summary of discussion session #1 - Willoughby 

Summary of discussion session #2 - Chan 

Summary of discussion session #3 - Hodur 

Break 

Field experiment status report - R. L. Elsberry 
Discussion of potential hypotheses 
Recurvature - G. J. Holland 
Recurvature - R. L. Elsberry 
Lunch 

Continue discussion of potential hypotheses 
Break 

Continue discussion of potential hypotheses 
Close 

30 January 1988 

Wrap-up session 
Coordination of write-ups 
Close 



37 



Appendix C 

Summaries or short abstracts of accomplishments during 
the first year of the ONR Research initiative have been 
submitted by Chi, Elsberry, Holland, Krishnamurti and 
Williams. DeMaria and Willoughby of the Hurricane Research 
Division (NOAA) submitted abstracts of their research that 
is relevant to this research initiative. 
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Heat, Moisture and Momentum Fluxes in the MPBL 
of Tropical Cyclones 

Joseph Chi, Ph.D. 

University of the District of Columbia 

Summary 

Under the ONR Grant, UDC has used state-of-the-art 
turbulence theory to simulate the marine atmosphere in 
tropical cyclones. Three main accomplishments are briefly 
described: 

1. Developed a mathematical model for the hurricane 
vortex flow over the sea surface. The model uses the 
second order turbulence closure equations. 

2. Developed computer program for calculating 
velocity, enthalpy, moisture, and turbulence energy 
distributions in the MPBL of the hurricane vortex flow 
over the sea surface. 

3. Developed a laboratory model of vortex flow over 
the water surface to compare calculations of turbulence 
theory with experimental data. 

Lists of publications supported by this research grant are 

also given. Finally, research needs in the future for 

conceptual advances to be made in the area of dynamics of 

marine atmospheres in tropical cyclones are presented. 

SCIENTIFIC ACCOMPLISHMENTS 

Second Order Turbulence Model for Hurricane Vortex MPBL 

A second order turbulence model has been developed for 
the marine planetary boundary layer of tropical cyclones. 

The model uses a cylindrical coordinate system. It permits 
three dimensional simulations of tropical cyclones that 
include influences of the centrifugal force, Coriolis force, 
air/water interface waviness and interface liquid droplet 
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entrainment. A complete set of conservation equations, 
Reynolds stress equations and turbulent flux equations for 
the three dimensional flow have been derived. Details of 
the mathematical model for the tropical cyclone are 
documented in a report entitled "A Turbulence Theory for 
Vortex Marine Planetary Boundary Layers." Some solutions 
for these equations have also been obtained. 

Computer Simulation of Vortex Boundary Layer Flows: 

A computer program has been developed to simulate the 
turbulent MPBL of tropical cyclones over the sea surface 
(Figure 1) . It includes the boundary layer heat, moisture 
and momentum transfer; turbulent energy production, 
dissipation and diffusion; strong influence of the 
gravitational, Coriolis, and centrifugal forces; influence 
of sea water entrainment in the vortex; and influence of 
air/water interface roughness. Figure 2 shows a set of 
typical results of distributions in the hurricane MPBL of 
the radial component of wind velocity, the tangential 
component of wind velocity, the moist static enthalpy, the 
total moisture ratio and the turbulent kinetic energy, 
respectively. 

Laboratory Experiments of Vortex Flows Over a Water Surface 

In order to increase confidence in the computer 
simulation, a laboratory model shown in Figure 3 for air 
vortex flows over the water surface has been designed and 
fabricated. The model's overall dimensions are 55 cm 
diameter by 90 cm high. The vortex is generated by forcing 
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air through 24 evenly spaced vanes of 15 cm high placed on a 
30.48 cm diameter pitch circle. The swirling air discharges 
from the model through a 3.81 cm diameter hole located at 
center of the model's top disk. The bottom of the swirling 
air is in contact with a pool of water that can be 
maintained at a constant depth and at desired temperatures. 
Figure 4 shows calculated distributions in the boundary 
layer over the water surface of the radial component of air 
velocity, the tangential component of air velocity, the 
static enthalpy, the moisture ratio and the turbulence 
kinetic energy, respectively. Measurements in the boundary 
layer have been planned using a laser doppler anemometer 
that is currently under development at UDC. Some 
preliminary experimental results of flows in boundary layers 
over solid surfaces (instead of water surface) have been 
obtained from measurements using the hot-film anemometer 
system. These data are plotted in Figure 5 for comparison 
with results of the computer simulation. Excellent 
agreement between theory and experiments can be seen in this 
figure . 
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FUTURE RESEARCH NEEDS 



It is suggested that both the UDC computer model for 
MPBL of tropical cyclones and the UDC laboratory model for 
the vortex boundary layer over the water surface should be 
expanded. For the computer model, the axisymmetrical flow 
assumption in the vortex simulation should be removed. 

Then, to increase calculation efficiency, the finite-element 
numerical scheme should be developed to replace the finite- 
difference calculation scheme. Instrumentation for the 
laboratory model should be improved to facilitate 
simultaneous measurement of air velocity, temperature, 
moisture ratio, and turbulent fluxes. For this purpose, 
optical instruments should first be developed to eliminate 
interference of instruments with the vortex flow. 
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FIG. 1 HURRICANE VORTEX MPBL SCHEMATIC 



FOR TURBULENCE MODELING 
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FIG. 2 COMPUTER SIMULATION RESULT FOR HURRICANE MPBL 





ENTHALPY (H-H q ) 





RADIAL INWARD VELOCITY (-U) flOISTURE (0-Q q ) 



i 



H 

E 

I 

G 

H 

T 

« 

K 

M 

# 




TANGENTIAL VELOCITY (v) 




0-0 o.o 



turbulence kinetic ENERGY (k) 



44 



AIA OUTLET 





FIG. 3 



1AB MODEL SCHEMATIC 



FIG. 4 COMPUTER SIMULATION RESULT UNDER LAB CONDITION 
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FIQ. 6 COMPUTER VS LAB DATA COMPARISON 



Predictability of Tropical Cyclone Motion 
Mark DeMaria 

Hurricane Research Division, AOML, NOAA 
4301 Rickenbacker Causeway 
Miami, FL 33149 



Abstract 

The dynamical method of the study of atmospheric predictability 
is applied to the tropical cyclone track forecasting problem. This 
method estimates the growth of forecast errors by comparing several 
numerical model simulations with slightly different initial 
conditions. The barotropic vorticity equation is initialized with 
vertically averaged winds from NMC analyses and is solved using the 
spectral method with spherical harmonic basis functions. The model 
uses triangular truncation at wavenumber 128 (about 100 km resolution) 
and 39 forecasts of Atlantic tropical cyclones which occurred between 
1979 to 1985 are considered. 

The effect of initial position errors is studied by running 
each forecast with the initial vortex center displaced from the 
observed initial position. These results show that initial position 
errors have little effect on the total forecast errors after about 
12 hours. This can be explained by geometric considerations and the 
fact that the distance between the displaced forecast track and the 
original forecast track remains small compared to the forecast error. 
The effect of initial position error becomes more significant if 
errors in the initial motion of the storm are introduced. 

The effect of inadequate data coverage is studied by comparing 
forecasts with various wavelengths filtered from the NMC analyses. 
Using this method, the lower limit to track forecast errors can be 
estimated for a given data resolution. These results show that 
inadequate data coverage introduces large errors in the forecast 
tracks. For Atlantic storms east of 70 W, the observed track forecast 
errors are close to the estimate of the lower limit of error, 
suggesting that little forecast improvement can be made without 
improved data coverage. However, for Atlantic storms west of 70 W, the 
observed errors are larger than the lower limit suggesting that 
further improvements can still be made. 
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Progress Report on Tropical Cyclone Motion Studies 

R. L. Elsberry 

Motivation from the DOD viewpoint for the Tropical 
Cyclone Motion Initiative included the Typhoon Abby (1983) 
case in the western North Pacific. Nearly all of the DOD 
assets were put on alert (unnecessarily in most cases as it 
turned out) because the operational forecast tracks were in 
error. The actual track of this very large and dangerous 
typhoon deviated significantly to the right of the "steering 
flow". One of the basic questions to be addressed was 
whether this anomalous track behavior could be attributed to 
the anomalous structure of the vortex. 

The Ph.D. dissertation of Mr. Michael Fiorino at the 
Naval Postgraduate School addressed this question. Fiorino 
used a numerical model to study a vortex structure that 
included stronger winds at large radii than are found in the 
normal typhoon. He also studied more intense inner-core 
wind effects to determine whether this aspect of Typhoon 
Abby was the cause of the anomalous track. Fiorino 's 
numerical tests demonstrate clearly that the strength of the 
winds in the outer regions (beyond 300 km) could account for 
the differences in track, whereas an increase in inner core 
intensity by 50% has a mininal effect on the track. The 
translation speed is almost linearly dependent on the wind 
strength in the outer (300-1000 km) radius of the vortex. 

This result has implications for the strategy of 
observing tropical cyclones and for proper initialization of 
the numerical track prediction models. First, observations 
are required at much larger distances from the storm than 
are generally available. Without aircraft reconnaissance, 
new remote-sensing techniques need to be developed and 
tested that will detect the differences in vortex structure. 
New Loran-based (or in the future, Global Positioning System 
based) dropwindsondes should be deployed in the outer 
regions of storms. Complete coverage over the required 
observational regions may not be achieved until the Laser- 
wind satellite system is available! Second, many of the 
present numerical track prediction models have a systematic 
deviation in the track that might be explained by Fiorino 's 
results. New observations to define accurately the actual 
outer wind structure need to be incorporated in these 
numerical models to improve the track predictions. For 
example, that part of the vortex motion (speed and 
direction) due to the so-called beta effect is closely 
related to the ventilation flow averaged over the inner 300 
km of the storm. 
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